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Abstract 

The narrowing of electron and ion wave packets in the process of photoionization is investigated, 
with the electron- ion recoil fully taken into account. Packet localization of this type is directly 
related to entanglement in the joint quantum state of electron and ion, and to Einstein-Podolsky- 
Rosen localization. Experimental observation of such packet-narrowing effects is suggested via 
coincidence registration by two detectors, with a fixed position of one and varying position of the 
other. A similar effect, typically with an enhanced degree of entanglement, is shown to occur in 
the case of photodissociation of molecules. 
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I. INTRODUCTION 



In the course of photoionization, a photoelectron is ejected and the ion recoils, being 
constrained both by the conservation of momentum and energy and by the condition of the 
original atom. In principle this initial condition includes all aspects of the atom's internal 
and center of mass states, but here we will focus for greatest clarity on an atom in its 
ground electronic state with a center of mass wave packet determined by whatever localizes 
the atom in the region of the photoionization. Although photoionization has been treated 
repeatedly in weak and strong fields (see, for example, Refs. Q and 3), the focus has 
been mainly on cross sections and the dynamics of the particles, and there has been little 
discussion of the nature of the joint quantum state of the breakup fragments. This quantum 
state is entangled and it is closely related to questions of fundamental interest because the 
fragmentation process is exactly the one used by Einstein, Podolsky and Rosen (EPR) |3|] 
to illustrate Einstein's position regarding limitations of quantum theory. We identify the 
amount of state entanglement by examining the relative localization of the wavepackets of 
the electron and ion, under the simplifying assumptions that the incident photon momentum 
and the post-breakup Coulomb interaction can be neglected. The different time regimes for 
entanglement are identified. 

Due to the finite mass of the particles, the wave function of the system changes drastically 
between the initial stage and the long-time limit. We find suitable measures of the entangle- 
ment of the electron and ion that are connected with their packet widths in position space, 
specifically the coincidence width and the single-particle width. In view of other treatments 
of breakup jj], these localization measures give an alternative view of entanglement and 
reveal new channels for achieving high degrees of entanglement. Our choice of measure also 
identifies entanglement "control parameters" for comparison with those that have been ad- 
vanced in previous studies of both photon-atom and photon-photon 6] wave functions, 
and through conditional localization it is related to formal photonic analogues 7 of the 
EPR discussion and bimolecular breakup as a route to matter wave entanglement q| . More 
importantly, these widths are experimentally measurable entities. Finally, we discuss similar 
issues for dissociation of a diatomic molecule and explain the most significant differences in 
the results. 
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II. PHOTOIONIZATION 



Let an atom, originally in its ground state, be photoionized by a light field 

S{t) = So sin(cut), (1) 

where huj > \Eq\ and Eq is the ground state energy. It should be noted that by using the 
dipole approximation and ignoring the term k - fin the argument in Eq. we are ignoring 
all recoil effects due to absorption of the photon momentum hk. This approximation is 
quite reasonable because there is another much stronger mechanism giving rise to recoil. 
In the process of photoionization an atomic electron acquires an energy ~ hu and hence 



a momentum ~ y/rnhuj, and the ion gets the same momentum (with the opposite sign), 
and this momentum is much larger than hk = huj/c. This is in contrast to the problems of 

n 

entanglement in spontaneous photon emission of excited atoms and Raman scattering p. 

To describe such a process with atomic recoil and with an initial wave-packet distribution 
of the atomic center of mass, let us begin from the Schrodinger equation for two particles - 
electron and ion - in the field. Traditionally, to separate variables in such an equation, we 
use the relative (rel) and center-of-mass (cm) position and momentum vectors 9| 

^ ^ ^ ^ mere + rriifi 

'rel ' e ' 1 5 ' cm s ,r ' , , 

Prel ) Pcm Pe ~^ Pit 

where rg and are the electron and ion position vectors, Pe = —ihd/dfe and Pi = —ihd/dfi 
are their momenta, and me and rrii are their masses, with M = rrie + mi. It is worth noting 
that the "mixed" coordinate- momentum variable pairs fVei and Pcm, as well as rem and Prei, 
each have zero commutator. For example, [fVebPcm] = 0. For this reason one can call them 
EPR pairs, recalling the famous discussion of Einstein, Podolsky and Rosen P| . 
The Schrodinger equation takes the form 

where /x = merrii/M is the reduced mass. Because we have made the dipole approximation, 
the variables r^ei and rem in Eq- © are separated, and its solution is a product of functions 
depending on these two variables separately: 

^(rrel, rem, t) = ^em(rem, t) X ^re^rrel, t) , (4) 



where the equations of motion of \l/cm('^cm, ^) and \I/rci(^rci, ^) are 
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dt 2/i rrci 

We note that the factorization shown in ^ is far from the same as factorization in the 
particle variables fe and fi. That is, the electron and ion are quantum entangled in the state 
given in Q. 

Let us assume that the initial atomic center-of-mass wave function is given by a Gaussian 
wave packet with width Ari2: 
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Then, as is well known P], the time-dependent solution of Eq. (0) has the form of a spreading 
wave packet such that 
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where Arcm(t) is the time-dependent width of the center-of-mass wave packet 
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and tip™^ is its spreading time, tsp™^ = 2M 



cm 



//i. At t > ti"^^ the width Arcm(t) grows 



linearly and the velocity of spreading equals to fip™'' = h/2MArim- 

Under the conditions of interest here, the solution of Eq. (jH} is only a little bit more 
complicated. The initial wave function of the relative motion is taken to be the hydrogen 
ground-state Is wave function 
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where -Rio('"roi) is the hydrogen radial wave function for the principal quantum number n = 1 
and angular momentum I = 0, and loo = l/"\/47r is the spherical function for I = rrii = 0. 



We assume a sufficiently high photon energy huj to ignore bound-bound transitions. Then 
the time-dependent wave function obeying Eq. (jH)) can be presented in the form 

^reKrrel, t) = Co{t) ijj,, + 6"^^* / dE Cs(t) ^Ep, (H) 

^0 

where ipEp is the field-free wave function of the continuous spectrum with / = 1, m; = 0: 

tpEpir.ei) = REi{rrci) yio(cos 6*1.61), (12) 
in which Rei{i^vci) is the radial wave function for energy E and angular momentum / = 1, 



Yio = ^y^/An cos^^rei, and 6'i.ei is the angle between the vectors Sq and f^ei- 

With multi-photon processes ignored, the equations of motion for the probability ampli- 
tudes Co{t) and CE{t) in the rotating- wave approximation, following directly from Eq. ©, 
are given by 

doE ■ £q 
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where dEo = (c^os)* are the bound-free dipole matrix elements of the atom, and we consider 
the case of a pulse with rectangular envelope, which means that the interaction is turned on 
suddenly at t = 0. 

With the help of adiabatic elimination of the continuum Eq. ()13ap can be reduced to 
a much simpler form: 

ihCo{t)-{Eo-tn-fj) Co(t) = 0, (14) 
in which the amplitude decay rate 7/ is half the Fermi-Golden Rule rate of ionization: 
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The solution satisfying the initial condition Co(0) = 1 is 



Co(t) =exp ( --Eot-^it ] . 



(15) 
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With this function substituted into the right-hand side of Eq. ()13b|) . the equation for CEif) 
can be easily solved to give, with the initial condition Ce{S^) = 0: 



CE{t) = - 
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At times t 'jj^ both Co{t) in Eq. and the first exponential term in Eq. p7j) vanish. 
As a result the wave function \l/i.ei describing relative motion takes the form 



^2 i>oo 

^ve\ir,eht) = — ^COs6'rol / dE REl{rj:el) 
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xexp(-^Et) ^ ^ ■ (18) 

We assume that the laser frequency io and hence the energy E ^ Eo + huj are high enough so 
that the radial function REi{fmi) is approximated by the well known high-energy field-free 
expression for the Coulomb radial wave function 
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where k = \^2fiE/h, Si is the Coulomb scattering phase for I = 1, and oq = h'^/fxe^ is the 
Bohr radius. 

When the photoelectrons have energy far above the continuum threshold, we have h^j -C 
E Eq + uj. In this way the lower limit of the integration over E in Eq. ()18j) can be replaced 
by — oo. The energy E is approximated by E^ = Eq + huj in all the pre-exponential factors 
except the denominator on the right-hand side of Eq. ()18|). Also both the scattering phase 
5i and logarithmic term in the argument of cosine in Eq. (fT^ are neglected, and the factor 
k in the product kr is expanded in powers oi E — E^, viz, k ^ k^ + {E — E^)/h,v, where 



cos 6'i.ei 
X exp 



/c* = v^2/ZE^//i and v = a/2-E'*/ jj, = fik^/ jj, is the velocity of the relative motion. Then the 
integral over E can then be evaluated by the residue method, giving 

-7/(t-^)]^^(^t-r.ei). (20) 

This equation describes a spherical wave packet in rj-ei with an angular modulation deter- 
mined by the factor cos6'i.ei, propagating in the direction of growing r^-d with velocity v, 
having a sharp edge at rrei = vt and an exponentially falling tail at rrei < vt. The radial 
width of the wave packet ()20j) is v/2'yj. 

III. EVOLUTION OF THE RELATIVE-MOTION WAVE FUNCTION 

Although we have assumed that the time t exceeds the total ionization time (t > 77^), 
Eq. (|^m) still describes the initial stage for the relative-motion wave packet evolution after 



ionization. In this sense is the initial width of the relative-motion wave packet |\I/reiP, 

which we denote with ^rfj = f/27/. This width can change later due to dispersion. To 
describe such a spreading effect, we can extend our series expansion of the function k{E) 
up to second order m E ~ E^: k ^ k^, + {E — E^)/hv — {E — E^)"^ /2hfiv^. This gives rise to 
an additional factor in the integral over the energy E: exp {—i E^rj.ci/2hfiv^}. To keep the 
possibility of integration by the residue method, we have to use the Fourier transformation 
of this factor 

E'r,,, 
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With this representation we first carry out the integration over E (by the residue method), 
and then the one over r. The result is 
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where Erf is the error function, and we have defined 

HTj-qI ^rel 
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spr — yu(Ar^^°[^)2//i is the spreading time of the relative-motion wave packet. As the 
value of |\I'reiP is concentrated around Tfei ~ vt, by putting Tfei ~ vt in the definition of the 
parameter (, we get ( = t/t'i^P . In this form the meaning of ( is obvious: it is the time after 
ionization measured in units of the spreading time of the relative-motion wave packet. In 
Figs, m^a) andmb) the function |\l/reiP is plotted in its dependence on p = (rrei — vt)/ Ar^^j 
at small and large C's respectively. 

In the small-spreading regime {( ^ 1) |\I'reiP returns to the form of Eq. (j^UI) . but with 
additional oscillations on the left wing and a slightly smoothed right wing as compared to 
the step function jump of Eq. (|^. In the large-spreading regime {( ^ 1) |\l/reiP takes a 
Lorentzian shape: 



where 

Arrei(t) = CArS = -^Ar^ = v^,,t (24) 

and Vspr = Ar^gj/tspr = ^/fJ'^^lci- Altogether, at small and large (, the time-dependent 
width of the relative-motion wave packet is given by 



Arrest) 



27/ 

Vsprt = = t > ti? (C > !)• 



(25) 



In spite of a difference between the Gaussian center-of-mass (jHl ) and relative-motion 
[ (PUjl - lj^ . (j2ni)] wave packets, their widths behave similarly: in dependence on t they start 
from the initial values Arcm and Ar^°j\ and at time t longer than the corresponding spread- 
ing time both Arrci(t) ()25j) and Arcm(i) © grow linearly. In both cases the spreading time 
is inversely proportional to the squared initial size and the velocity of spreading is inversely 
proportional to the initial size to the first power. The only qualitative difference concerns 
the mass of an objet: the total mass M of the center-of-mass wave function is substituted 
by the reduced mass /i in the case of the relative- motion wave packet. 

The relation between the center-of-mass and relative-motion wave packet widths can 
change with time due to different spreading velocities of these wave packets. This makes the 
time evolution of the electron-ion wave function rather complicated, and this problem will 
be discussed separately in Section IVl 



IV. LOCALIZATION OF THE ELECTRON-ION WAVE PACKET AND ENTAN- 
GLEMENT 

Before going further into details of time evolution, let us discuss the entanglement effect. 
As we will show, in cases of initially localized pairs of particles as in photoionization, and 
where significant further interaction is absent, entanglement can be evaluated by carrying 
out a series of localization measurements. This has a close analog in earlier studies of 
.po^aneous photo, emission w.h a.o. .ecoi. Q B Q a. we., a. in the ■neas.e.ent- 
induced localization and entanglement discussed recently in a very different context by Rau, 
et al. 
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We will proceed by determining the dependence of entanglement on the ratio of widths 

- (26) 

where the time t is taken as a parameter. Note that in our treatment r]{t) is constrained 
only by momentum and energy conservation (e.g., we ignore final-state electron- ion Coulomb 
effects). Here Arcm(^) is of kinematic origin whereas Arj.ci{t) is due to the dynamics of the 
ionization process. We will see that r]{t) acts as the sole control parameter for entanglement 
of the two-particle system. 

In accord with Eq. (jH), the product of the wave functions Eqs. (jHl) and (pT|l determines 
the total wave function of the ion-electron system. It should now be considered as a function 
of ion and electron position vectors 

^(fe, fi, t) = \^ , t X ^rel(re - ^i, t) , (27) 



M 

showing that both \1/ and its squared absolute value are not factorable in the individual 
particle coordinates fg and fi. Such non- factorization defines quantum particle entanglement 
of electron and ion. 

Now we focus on measurements appropriate for seeing entanglement. We need to dis- 
tinguish coincidence and non-coincidence (single-particle) measurements, which have their 
theoretical counterparts in conditional and non-conditional probability distributions. For an 
example of a single-particle measurement, the electron probability distribution is measured 
regardless of the ion position (or vice versa). In contrast, a coincidence measurement assumes 
that a distribution of electron positions is registered while the ion detection position is kept 
at a given (constant) location (or vice versa). The difference between the results of coinci- 
dence and single-particle schemes of measurements is illustrated by Fig. |21 In this picture, 
in one dimension, we shade the region in which the joint probability density \'${fe,fi,t)\'^ 
is significant. In the left plot the sharp leading edge of the theta function in Eq. (j^Uj) is 
apparent, with its long exponential tail, and one also sees the more abrupt Gaussian cut-off 
on the sides. A purely schematic view of the same thing is shown in the right plot, where 
artificially sharp dashed-line borders are introduced, and supposed to be determined by the 
localization zones of the relative-motion and center-of-mass wave functions. 

Consider first an examination of the electron wave packet by the coincidence-scheme 
method, for a given ion coordinate Xi = const. The normalized measure of its width, 



Aa;e/Ari.ci(t), will be given by the distance between the points marked a and h. In contrast, 
the single-particle width takes into account the contributions from all possible different Xj's. 
Thus a suitable measure of the single-particle width of the electron wave packet is given by 
the distance cd. It is obvious that the electron packet is relatively highly localized when 
cd ^ ah. Correspondingly, a horizontal line through the shaded region would provide a nor- 
malized measure of Axj, etc. From this sketch we formulate two conditions simultaneously 
necessary for entanglement to be large: a high aspect ratio of the shaded area and a nearly 
diagonal angle between the dashed lines restricting the wave packet localization zones and 
the coordinate axes Xe and Xi. The high aspect ratio condition means that one of the two 
wave packets ("cm" or "rel") is much wider than the other one. 

Now we note that this relatively great localization condition is the same as a high en- 
tanglement condition. This becomes obvious by considering the two-particle wave packet 
|\l/(re, Tj, t)p as an information container 2. Then the bibliography gets the added element: 

Entanglement means that knowledge of one of the particles imparts information about 
the other, whereas non-entangled particles provide no information about each other. The 
greatest cross-specification of joint information by an entangled packet occurs when knowl- 
edge of one particle automatically corresponds to precise information about the other, i.e., 
knowledge of position Xi strongly localizes the region where Xe can be found. Reflection 
shows that a thin diagonal packet in Xi-Xe space achieves this, and normalized relative in- 
formation gain is well expressed by the ratio of single to coincidence width. Note that left 
or right inclination of the diagonal is immaterial. In contrast, in wave packet language the 
non-entangled condition (information about Xi gives no information about Xe) is equivalent 
to a factored wave packet: \E'(xj,Xe) 4>i{xi) 02(a^e)- One sees that a sketch corresponding 
to Fig. 121 but for independent particles, has dashed lines that are horizontal and vertical, in 
which case all Xj's predict exactly the same Ax^, and vice versa. 

Mathematically, single-particle probability densities are given by Eq. (j27jl integrated ei- 
ther over fi or fg: 

Pe{re,t) = j dfi |^(r;,r;,t)|2, (28) 

or 

Pi{ri,t) = j dr, \^{f,,ri,t)\^. (29) 
Such distributions reveal no entanglement effects because all the information about the 



10 



position of one of the particles is lost completely when the two-particle probability density 
is integrated over or fe. However, Pe and Pi do serve a normalization role, as we explain 
later. 

Let ^r^f^ and ^rf'' be the widths of the single-particle electron and ion wave packets, 
where, for example, |Are''^p = (jr'eP) — K^e)Pi with 

(^e) = j dfe fe Pe{fe,t) 

dfedfi fe |\I^(r;,fj,t)|^ 



and 



TTl ■ 

{rcm) + J^{rrel) (30) 



(|fe|') = j dfe rl Pe{fe,t) 



dfedfi rl |^(r;,r;,t)p 



rrii 



2 



— ( ^cm + j^^i^cl ) 

= (Ir-^l^) + 2^(r-;^)(r-;ei) + ^drVeil'). (31) 

Note that we have used the relation fe = fcm + j^r^ei and changed the integration variables to 
the center-of-mass and relative coordinates. Then Eqs. ()3()|1 and (p?T|) yield the single-particle 
measures: 

and, similarly. 



Note that the relative- motion wave packet width Arrei(t) plays the role of a natural normal- 
ization factor for both single-particle and coincidence-scheme (see below) electron and ion 
wave packet widths. Divided by Ari.ei(t) these widths become dimensionless, and they are 
denoted dr^tj and Sr^^j. 

In the coincidence scheme of measurements the overall width of the distribution ()27j) with 
respect to fe at a fixed fi is given by the smaller of {M/me)Arcm(t) and Arrest), which is 
well-represented by a simple formula for the coincidence measures: 

Sr'f> . ^ « , (34) 
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The expression in Eq. ()34p is a better approximation the closer we are to one of the extreme 
cases (M/me)Arcm(i) ^ Ari.oi(t) or ^ Ar^ci{t)- As shown later in this section, these limits 
correspond to the high entanglement regimes of main interest, and we do not need to bother 
too much about the details of the intermediate region. Similarly, at a given re, the widths 
of |\l/cmP and |\l/reiP with respect to rj are correspondingly (M/mj)Arcm(^) and Arrest). 
Therefore the overall width of at a given is the smaller of (M/mj) Arcm(t) and 
Arrei(t), corresponding to the formula 

fr<=' = ~ "W ,35, 



.MJ 

It is clear now how rjit) of Eq. (j^Bj) serves as a "control parameter" for both coincidence 
widths. Plots of 6rf^^ and 6r^'^^ as a function of r/, are shown in Fig. Efa) whereas graphs of 
6rf^^ and 5rj'^'' are shown in Fig. IH^b). Note that we use in these graphs an artificial value of 
the electron to ion mass ratio nie/mi = 0.1 so as to show more clearly the difference between 
the two curves for r/ < 1. However, all the qualitative conclusions from these pictures remain 
the same for a more realistic value of this ratio m^/mi ~ 10~^. One of these conclusions is 
that we always have dr^f'^ > Sr^^}. 

The ratios of single-to-coincidence electron and ion wave packet widths, dr^^j to Sr^^j, can 
be considered as a measure of entanglement, as remarked at the beginning of this section. 
As we will show elsewhere they are essentially identical to the corresponding Schmidt 
number discussed in earlier discussions of photon-atom entanglement jj,!^. These ratios 
can be referred to as the electron and ion entanglement parameters in the form: 



= ^ and = (36) 

Entanglement is large if -Rg ^ 1 and/or Ri ^ 1. If i?e ~ 1 and i?j ~ 1, there is little or no 
entanglement at all. 

By using Eqs. (jH^ - (jHSl), we can find a useful approximate form of R^ and Ri : 



which is plotted in Fig. ^ Note that here the width ratios for the electron and the ion are 
the same, i.e.. Re = Ri- This is actually true only when the widths of \l/cm and \l/rei are very 
different from each other, or equivalently r/ ^ 1 or r/ ^ 1. Even though R^ and Ri may not 
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be exactly the same in the zone nie/M < t] < 1, they both have values around unity, which 
corresponds to the relatively less interesting small entanglement regime. Thus we designate 
the two of them together by R without a subscript. The asymptotic behaviors of R in three 
different regions of ri{t) are particularly noteworthy: 

Region 1, V«j^«j^: 

mem, Ar,ci(t) f fi \ 1 , . 



M2 Ar,^(t) \MJr]{t) 
Region 2, ^ « « ^ : ~ 1; (39) 

Regions, ^ « ^ « r/ : Rr^r^it). (40) 
Note that the minimal value of the entanglement parameter pUj) is equal to one, -Rmin = 1, 



and it is achieved at ri{t) = a/ yu/M. 

V. TIME EVOLUTION OF PACKET WIDTHS AND ENTANGLEMENT PA- 
RAMETER 

In Figs. El and 1^ both the electron/ion wave packet widths and the entanglement param- 
eters are shown, for fixed t, in their dependence on the control parameter defined earlier: 

However, we can use the same pictures to show the time evolution of the widths Are(t) and 
Arj(t) and the entanglement parameter R{t) defined in (jSZI)- To do this, we have to learn 
how ?7(t) changes with time. 

The two typical cases of significantly different behavior are illustrated in Figs. El and |B1 
Parts (a) and (b) of these Figures show the time dependence of the widths Arcm(^) and 
Arrci(t) themselves and of their ratio, which equals ri{t). A key feature of ri{t) is its strong 
dependence on the initial sizes of the center-of-mass and relative-motion wave packets, AriS 
and Ar^gj, or in other words, on the initial value of the control parameter 77(0) = rjQ. 
Depending on its initial value, rj is either rising as shown in Fig.jSl^a) or falling as in Fig.jHfb). 
The border between these two regimes is given by rjo = rj^, where 
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If rjo < rj^,, the center-of-mass wave packet spreads faster and eventually becomes wider 
than the relative-motion wave packet, though initially Arcm ^ ^'"^d • this case the control 
parameter ?7(t) is a monotonically growing function of t (see Fig. Efb)). On the contrary, 
if rjQ > 1]^, the center-of-mass wave packet spreads slower than the relative-motion wave 
packet. Though initially the center-of-mass packet can be either narrower or wider than the 
relative- motion packet, at very large t the relative- motion packet becomes wider than the 
center-of mass packet. This gives rise to a falling function //(t) shown in Fig. Efb). In both 
cases (?7o < rj^, and rjQ > rj^,) the ranges of variation of the parameter 77 (t) are finite. At very 
long times 77 (t) has the asymptotic value 

" ^ - Jih («' 

which follows directly from the definition of 1] (j41|) and Eqs. Q and (j25p for the widths 
Arcm(i) and Arj-dit). In the case rjo = rj^ the parameter rj{t) does not depend on time at 
all: Tjit) = const (= rj^ = r]o = rjoo). 

By finding the evolution regimes for the control parameter //(t), we can also make defi- 
nite and interesting conclusions on the evolution of the entanglement parameter R{t) ()37p. 
Directly from Eq. (jHTj) one can easily see that for an arbitrary value of rj, the entanglement 
parameter obeys the relation 

R(-^)^Riv)- (44) 



The initial and final values of r] are connected with each other exactly by the same sub- 
stitution as used in Eq. we see that the initial and final values of the entanglement 
parameter must be equal: 

Rq = R^ = R{t 00). (45) 

For the entanglement parameter given by Eq. (j37p this equality is valid identically for all 
values of rjQ. If rjo is located in one of the high-entanglement regions of Fig. |3J rjo <^ fi/M 
as in (jHHj) . or ?7o ^ 1 as in (HOJ, the final value of r] is in the opposite of these two high- 
entanglement regions, rj^ \ 01 rj^ <^ ii / M . 

Thus we see that the time-dependent entanglement parameter R{t) starts from a large 
value i?0) falls to i? ~ 1 and then grows again to the same value from which it started. 
Physically such an evolution means that initially one of the wave packets is much wider 
than the other one, and for this reason the electron-ion entanglement is large. Then, as the 



14 



narrower wave packet spreads faster, they become approximately of the same width, and 
this corresponds to a small entanglement. Finally, when the initially narrower but faster 
spreading wave packet outstrips the initially wider but slower spreading one, the relation 
between their widths reverses, and this returns us to the case of a large entanglement. 

The difference between the cases ?7o ^ /i/M and r^o ^ 1 concerns only the direction of 
evolution, correspondingly, to the right or to the left in the r^-axis in Fig.|3J If the initial value 
of the parameter rj is located in the small-entanglement region (jlOJ, fi/M < rj < 1, all the 
conclusions about the direction of evolution and about the relation between the initial and 
final values of the entanglement parameter remain valid. However, in this case, at all times 
t the entanglement parameter remains on the order of one. If t]q = rj^, the entanglement 
parameter does not change at all, and R{t) = 1. 

VI. EXPERIMENTAL CONSIDERATIONS 

The discussion as given so far doesn't treat some elements that will come into play in 
experimental tests. In order to bring them into focus briefly, we show in Fig. [7| what can 
be called experimentally realistic zones. We have plotted the region where the relative 
probability distribution is non-zero. It has a three-dimensional aspect that we do not need 
to show because it is axially symmetric about the polarization axis of the ionizing light 
beam, taken as the vertical axis here. It is not spherically symmetric because of the dipole 
character of photoionization (i.e., the factor (Ie^q ■ Sq in Eq. ()20j)). The new-moon shaded 
areas indicate the regions where the relative-motion wave function I^E'reil^e)!^ is relatively 
large. 

Since the time evolution of the relative wave function is strictly limited by the step 
function 6{vt — rrei), we will here consider the ion position to define an origin of polar 
coordinates (r^ = 0), in which case a circle of radius = vt limits the range of the electron 
coordinate at time t. The relative coordinate probability distribution I^E'reK^e)!^ is of course 
not uniform inside this circle, so we have drawn the boundary on which |\l/reiP equals ^ of its 
maximum value. This creates two sectors with "new-moon" shape where there is the highest 
probability to find the electron, given that the ion is at the origin of the circle, and taking 
only |\E'reiP into account. However, the probable position of the electron is also influenced by 
|\E'cm(^e)P at Tj = 0. In the figure the black dot regions show the range of electron positions 
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given by |\l/cm(^e)P at different positions of the center of mass along the rg-rj hne. 

The figure has many variations, and the sizes of the new-moon and Are zones change in 
time, as our formulas indicate. The overall shapes will remain the same, and a generic high- 
entanglement experiment will be one that ensures overlap between a small black-dot region 
and a high-probability portion of one of the new- moon regions (near the circle boundary). 
Higher count rate, although reduced entanglement, will be associated with increased size of 
the black dot region. 

VII. PHOTODISSOCIATION 

It is easy to see that very similar results will arise in a treatment of photodissociation of 
molecules. Here we remark briefly on some of the differences. Let us assume that we consider 
a diatomic molecule undergoing dissociation. There will be a relevant dissociation rate 7^), 
which can be substituted for the 77- governing ionization, and just as for the atom there 
will be an initial localization of the molecular center of mass. Then the main differences to 
the ionization example arise because the mass ratio of the fragments is much closer to 1. 
Compared to the case of photoionization, where rrie <^ rrii, the masses Mi and M2 of the 
photodissociation fragments obey Mi ~ M2. Given this, the relative-motion velocity after 
dissociation v ~ ^jTiuj/ ^ is significantly smaller than in the case of ionization (in atomic units 
fmoi ~ a/^^/M). The main difference between photoionization and photodissociation results 
concerns the region < iri{t) < 1 of intermediate values of Arcm(^) in Fig. |3] where i? ~ 1. 
For Ml = M2 = \M this region degenerates into a single point rj{t) = ^ = 77* ()42|) . The 
entanglement coefficient R is large both at rjit) < t]^ and rjit) > r]^. To show more clearly 
the difference between photoionization and photodissociation we plot in the right picture of 
Fig. IHlboth molecular and atomic entanglement coefficients in their dependence on ln(?7), 
with the electron to ion mass ration taking a realistic value rrie/mi = 10~^. This picture 
shows that if in the case of photoionization there is a rather large region of intermediate 
values of rj where entanglement is small, i? ~ 1, in the case of photodissociation of a molecule 
the entanglement parameter is large practically at any rj except one point rj = rj^,. 

In dependence on time t the control parameter ri{t) changes in a way similar to that 
described above for photoionization: ri{t) grows if initially it is small {rjo < 77*) and falls if 
large (r/o > V*)- The final value of the control parameter r/oo is related to rjo by Eq. (jlBI), 
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which takes the form t]oo = 1/{4:T]o). As shown previously, the initial and final values of the 
time-dependent entanglement parameter R{t) are equal to each other, Rq = i?oo- At rjo = r/* 
both the control parameter ri{t) and the entanglement parameter R{t) do not change with a 
varying time t, ri{t) = rjQ and R{t) = 1. This is the only case when there is no entanglement 
at any time t. In all other cases (t^q 7^ V*) ^he entanglement parameter is large initially, 
reaches R = 1 one at such t that gives rjit) = rj^,, and then grows again until it reaches its 
initial value Rq- 



VIII. CONCLUSION 



We have evaluated the space-time behavior of the joint quantum state of an ion and 
electron following photoionization. Neglect of the incident photon momentum and of the 
final state Coulomb interaction means that the evolution of the state, and thus of the entan- 
glement between the two particles, is constrained only by free-particle two-body momentum 
and energy conservation. This evolution provides an exactly calculable illustration of the 
situation involving massive particles sketched in the famous paper of Einstein, Podolsky, and 
Rosen P]. We have obtained expressions for the entanglement-induced wave packet narrow- 
ing that occurs, and have indicated how entanglement can be identified and determined 
quantitatively. To do this we introduced R, the ratio between the entanglement-free wave 
packet width and the coincidence wave packet width. This is essentially the degree of entan- 
glement. We gave expressions for R in terms of ionization rate and packet spreading velocity, 
which are of course themselves determined by underlying parameters such as atomic bound- 
free dipole moments, relative electron and ion masses, ionizing field strength, etc. It was 
shown that R depends in a simple way on the basic control parameter rj = Arcm(i)/Ari.ei(t), 
and can be much larger than unity in two limits, when rj ^ 1 and also ^ 1. The same 
formalism can be applied equally well to photodissociation of a diatomic molecule. For 
realistic physical values of the relevant parameters, in a typical example of atomic photoion- 
ization, R is not very large because of the extreme discrepancy between rrii and me, but for 
photodissociation of a diatomic molecule, where the fragment masses can be approximately 
equal, R can be substantially increased. 
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FIG. 2: Views of the one-dimensional equivalent of in Eq. (|27jl . Here we illustrate the relation 
between entanglement and distribution of particle positions. High correlation between Xe and Xi 
can be attained only under the condition that the ranges of available Xe and Xj are large compared to 
the variation range of one of them at a fixed value of the other variable. We have used me/rrii = 0.2, 
r] = 0.5, and 'jit = 4 for illustration. The dashed lines drawn in the figure tell approximately the 
region of localization of the wave function. 
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FIG. 3: Electron (a) and ion (b) wave-packet widths in the schemes of single-particle and coinci- 
dence measurements with rrie/mi = 0.1. 
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FIG. 4: Plot of the entanglement parameter i? as a function of r/(t) with m^jmi = 0.1; ry* is 
the stability point (|42j) at which 7]{t) = const = y/Ji/M. The insets give the corresponding plots 
of the one-dimensional analog |^'(xe, a^j, t)P from Fig. |5J The axes of the three insets have been 
rescaled so as to show the details more clearly. The large entanglement regions are clearly seen to 
correspond to a large aspect ratio of the shadowed areas. 
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FIG. 5: Part (a) shows the time-dependent widths of the center-of-mass and relative-motion wave 
packets (in units of Ar^^j), and part (b) shows the control parameter ri{t). We have taken r/o = 0.05 
and rUe/M = 0.1. 
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FIG. 6: The same as in Fig. |21but with r/o = 0.5. 
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FIG. 7: The new- moon shaded areas indicate the regions where the function |^'roi(^e)P is relatively 
large, relatively large, i.e., at least as large as one-third the maximum value. Black dots indicate 
regions where |^'cm(^e)P 7^ 0, at a given rj. Three shown experimental situations correspond to 
different locations of the electron detector, relative to the ion position, which defines the origin. A 
smaller size of such a black dot inside the shaded area corresponds to a higher level of entanglement. 
In the case when the black dot is located far outside of the shaded area, there is no overlapping 
between the center-of-mass and relative-motion wave functions and the total two-particle wave 
function equals zero, ^' = 0. 




FIG. 8: Entanglement parameter for two dissociating molecular fragments with Mi = M2 (left) and 
the same dissociation curve plotted vs. ln(?7) on the right, where the corresponding photoionization 
curve is included for comparison, with its very different mass ratio, rrii = lO^mg. 
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